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^ ■ Abstract 

The conditions for the appearance of a sharp laser transition are formulated in terms of a scaling 
limit, involving vanishing cavity loss and light-matter coupling, k — V 0, g — v 0, such that g 2 /k 
stays finite. It is shown analytically that in this asymptotic parameter domain, and for pump rates 
above the threshold value, the photon output becomes large in a sense that is specified, and the 
photon statistics becomes strictly Poissonian. Numerical examples for the case of a two-level and 
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a three- level emitter are presented and discussed in relation to the analytic result. 
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I. INTRODUCTION 

An early preoccupation of laser theory was the analogy between the onset of lasing and 
phase transitions [lH3j. Within the cavity QED models the problem became increasingly 
accessible to accurate numerical treatments, and evidence that an abrupt change of regime 
takes place has accumulated. In parallel, approximate analytic considerations also argued 
in the same sense. 

A particularly illuminating case is the so-called random injection model of Scully and 
Lamb [4], which has the advantage of addressing directly the cavity mode statistics. The 
model was extensively studied 5j and is by now textbook materially, |7|. In emitter-plus- 
mode models the situation is more complicated because the photonic state information has 
first to be extracted, by eliminating the emitter degrees of freedom. This procedure can 
be carried out|8| in the case of single emitters with few (usually only two) states. Then 



information about the photon statistics can 
simplifying approximations, analytically too 



je obtained, either from the numerics or, using 
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In this context the seminal paper by Rice and Carmichael, Ref. [ill, has drawn the 
attention on the necessity of a limiting procedure for obtaining a sharp transition, with a 
precisely defined threshold. The analogy with the thermodynamic limit in the theory of 
phase transition was invoked. It was argued that for the lasing transition the limit involves 
the /3-factor going to 0, together with the cavity loss k so that the ratio /3/k remains finite. 
In this limit, and for a pump rate exceeding a threshold value, the number of photons N 
becomes infinite, generating "an explosion of stimulated emision" 111], and the appropriate 
object of study is the rescaled value (3N, which stays finite. 

In the present paper we show that a similar scaling limit is needed for an abrupt onset of 
lasing, but we use in our formulation the more ubiquitous Jaynes-Cummings (JC) coupling 
constant g, instead of the /3-factor. The latter is proportional to g 2 and indeed our scaling 
procedure implies k — y 0, g — y so that g 2 /n is finite. In this scaling limit and above 
the threshold one obtains N — y oo but such that kN remains finite. Moreover we are able 
to prove not only that the transition becomes abrupt, but also that the photon statistics 
above the threshold turns exactly Poissonian. The proof is analytic and does not rely on 
approximations. Also, analytic expression for relevant data, like the threshold pump rate, 
level occupancies and photon output in the lasing regime, are obtained. Numerical results 



are shown as illustration of the statements. 

II. THE MODEL AND STATEMENT OF THE RESULT 

Single emitter lasers are commonly described as embedded JC systems. In other words, 
two emitter states, either quantum dot or atomic configurations, interact with the cavity 
mode via the JC Hamiltonian 

H JC = gtf\2)(l\+gb\l)(2\ , (1) 

in the presence of, possibly, other states. Here b, tf are the photonic operators and the emitter 
states are denoted by \i). In particular |1) and |2) specify the upper and the lower laser 
state, respectively. We assume that the cavity mode is resonant with the laser transition. 
Dissipation effects are included in the master equation for the density operator p (in the 
interaction picture and with h = 1) 

l p =-i[H JC ,p]+£(p) (2) 

by the Lindblad terms 

£(P) = 2 [ 2b P b] ~ b ^ b P ~ P^b] +5Z ~2 2a v pa l ~~ a h a ijP ~ P a ij a iJ ' a ii = I*) 01 • ( 3 ) 

(Hi) 

A central role here is played by the first term, describing the cavity losses at a rate k, while 
the second term summarizes transition processes from states \j) to \i) at rates 7^-. Lowering 
a,ij operators and their hermitian conjugates <tJ- give rise to Lindblad terms accountin g fo r 



m, 



relaxation, but raising terms are considered as well, to simulate incoherent pumping 
and the corresponding rate will then be denoted by P. 

The master equation is solved in time until a steady-state solution is reached, from which 
data concerning level occupancies and photon statistics can be extracted, as function of the 
pumping rate and other parameters. With this information at hand one can address the 
problem of laser transition: how to identify it and what are the conditions for its appearance. 

The answer to the first question is that the lasing regime is defined by accumulation of a 
large number of photons in the cavity, the statistics of their number n obeying a Poissonian 
law 

Pn,n = ^2 Pn,n = ~ 6 ~ X ' for &11 H " ( 4 ) 



Equivalently, the Poisson statistics amounts to the requirement that the normal-ordered 
expectation values p n = (b^ n b n ^ depend exponentially on n, p n = X n , or that the zero time- 
delay n-th order correlation functions g^ = p n /Pi all become equal to 1. The problem is 
that, on the one hand, the large number condition is imprecise (how large is large?) and, 
on the other hand, no matter which of these criteria for Poisson statistics one chooses to 
apply, one has to check an infinite set of equalities. This is practically impossible, either 
experimentally or numerically. This is why it is often encountered in the literature that one 
limits oneself to simpler lasing criteria, like N > 1 and g^ = 1. 

The aim of the present paper is to formulate and prove the conditions under which 
strict Poissonian statistics is generated and at the same time to specify in what sense the 
photon output becomes large. Essentially these conditions involve the limit k — > and 
simultaneously g — > 0, but with the JC coupling parameter scaling like a/k so that the ratio 
g 2 J ' k remains finite. The necessity of a certain limiting procedure for obtaining a well-defined 
transition to a purely Poissonian statistics was recognized long ago|ll| and is analogous to 
the thermodynamic limit in the theory of phase transitions. The precise formulation of our 
statement runs as follows: 

a. Rescale k as e k and g as \feg. Then, in the limit e — > and for the pump rate above 
a threshold value, the average photon number N = p\ tends to infinity in such a way that 
the product sN remains finite. 

b. In the same limit and above the threshold, the rescaled expectation values p n = e n p n 
remain finite for all n, with their limit values obeying an exponential law p n = p™ . Equiv- 
alently, all the correlation functions become g( n > = 1. Below the threshold the increase 
in photon production is not sufficient, leading to p n = e n p n — > for all n > 1, and the 
correlation functions are in general different from 1. 

c. In this scaling limit the transition between the two regimes becomes sharp, with a 
well-defined threshold point. 

Several comments are in order. The model parameters, g and /t,that are rescaled corre- 
spond to the photon source and sink, respectively. They are present in all laser models and 
therefore the formulation of the scaling limit in these terms makes sense in various situations. 
The limit of small cavity losses (good cavity) pleads in favour of photon accumulation, but 
with the simultaneous diminishing of the production rate, proportional to g 2 , their overall 
increase is not a foregone conclusion. The statement (a) above spells out what is meant by 



a large photon number, namely that it should scale like 1/ ' k to compensate for the reduction 
of the rate of escape from the cavity. Thus, even for a cavity quality factor Q close to infinity 
there is still light coming out from the device. 

In situations as those described by incoherent excitation, at high pump rates the phe- 
nomenon of self-quenching[16] might occur. Beside producing population inversion, the 
pump can destroy coherence between the laser levels and inhibit the transition. This plays 
against lasing and therefore one may encounter a double transition, one at the onset of 
lasing and the other when lasing becomes quenched. In such cases the lasing regime takes 
place in a given interval of the excitation rates, limited below by the threshold value and 
above by the self-quenching P t h r < P < P sq . The endpoints of this interval depend on the 
model parameters. Accordingly, it may occur that the interval shrinks to zero and then no 
transition takes place. Such situations will also be discussed below. 

In what follows we will first bring numerical evidence in favor of the scaling limit result. 
We illustrate the situation with calculations performed on a two-level and on a three-level 
model. In both cases the scaling limit tendencies are quite clear. Still, numerical statements 
do not amount to a proof, which can only rely on an analytic argument. We are able to 
formulate an analytic derivation of the scaling limit result in the two- level case (see Sec. IIIII ) . 
Also, an analytic proof is available|17] for the random injection model[4|, [5|, which does not 
belong to the class defined by Eqs. (T5J), ([3]). We believe that all these results speak in favor 
of a wider generality of the scaling limit statement. 

A. Numerical results 

Steady-state results are obtained from the long-time limit of the master evolution. 
Particularly indicative of a transition is the behavior of the population inversion w = 
( |1) (1| — 1 2) (2 1 ) as a function of the pump. Depending on the parameters, one clearly 
detects two types of behavior |l2|, [14| illustrated in FigJjTa) 18]. The plots correspond to 
a two-level emitter in which the pumping is described by the raising |2) — > |1) Lindblad 
operator with the rate P = 712, while for the rate of loss to non- lasing modes we use the 
notation 721 = 7. For one set of parameters the curve is strictly concave, while for the other 
there appears an almost perfectly linear shortcut [9l. 1 171. 1 191] separating two concave regions. 



This clearly suggests that in the latter case an abrupt change of regime is taking place, in 
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FIG. 1. (a) Population inversion w in linear and (inset) semilogarithmic plot, and (b) second 
order correlation function g^ 2 ' for a two-level model. The parameters are, red (solid) line: 7 = 0.02, 
g = 0.1 and k = 0.01, blue (dotted) line: 7 = 0.01, g = 0.01 and k = 0.02. 



a given P-interval. It is also obvious that the appearance of the transition is conditioned 
by certain parameter values. By examining FigJT](b), we see that this linear segment (which 
in a semilogarithmic representation appears as an additional convex region, see inset of 
FigflK), corresponds to g^ 2 ' being very close to unity, which is characteristic for coherent 
light. Therefore it is natural to assume that we are in the presence of the lasing regime. 
Moreover, the fact that for large P-values the linear behavior disappears is consistent with 
the inhibition of lasing by self-quenching. 

Zooming in on the leftmost point of the linear interval, as in Figj2j it is seen that the 
transition becomes more and more abrupt as the scaling parameter becomes smaller, in 
accordance with the scaling limit statement. This is seen both in the panel (a) of Figj2j 
which refers to the two- level model discussed in Fig {H and in the panel (b) , which shows the 
result for a three-level emitter. In this latter case the pump is raising the system from the 
lower laser state |2) to a third state |3) (P = 732), wherefrom it relaxes to the upper laser 
state |1) with the rate 7' = 713. As before 721 = 7. 

The same tendency is seen in plots of g^ 2 \ which becomes closer and closer to the coherent 
light value g^ = 1, as the scaling parameter becomes smaller, both for the two- and for the 
three-level case, as illustrated in FigJU 

Finally, in FigJUwe show numerical results for the photon output. It is seen that in the 
lasing interval the rescaled photon numbers N = eN have practically reached their limit 
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FIG. 2. Scaling parameter dependence of the population inversion w for (a) the two-level emitter 
with parameters: 7 = 0.02, g = -y/eO.l and k = eO.Ol and (b) the three-level emitter with 
parameters: 7 = 0.02, 7' = 0.05, g = y/sQ.l and k = eO.Ol. Thin solid lines correspond to the 
analytic result Eq. (fT7|) . 
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FIG. 3. Scaling parameter dependence of g^ 2 ' for the same parameters as in Figj2j 



values, given by Eqs. ( 1T8|) and (TT9|) . also plotted in the figure. This means that indeed, iV 
grows like 1/e, in accordance with the scaling statement. 

The agreement seen in this section between the numerical data and the scaling limit 
results can hardly be accidental. Values of g^ close to unity suggest that in the interval of 
intermediate pump strengths the system operates in the lasing regime, and the large photon 
output speaks in favor of this supposition too. Nevertheless, it is not at all clear yet how 
the linear dependence of the population inversion on the excitation is in any way linked 
to lasing. The analytic arguments of the next section will prove that, indeed, the two are 
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FIG. 4. Rescaled photon output for the two- and three-level emitter cases. The parameters are 
the same as in Figj2j Thin solid lines correspond to Eq. (|18p in panel (a), and to Eq. (|19p in panel 

(b). 

related and appear simultaneously. 

III. THE TWO LEVEL LASER 

A special feature of Eq. fl2]) is the fact that it provides a system of closed equations for 
a subclass of relevant density matrix elements. In this category, the elements which are 
diagonal in the emitter states \i) are also diagonal in the photon number n, /?Jf n , and the 
only off-diagonal elements are of the form p n ' n+1 and their complex conjugates p n ' +1 n . This is 
a consequence of fact that the JC Hamiltonian conserves the excitation number |1) (1| + tfb. 

We consider here the case of an emitter consisting of only the two laser states |1) and 
|2). The Lindblad terms describe, beside the cavity losses, the spontaneous emission into 
non-lasing modes and the pumping, with the rates 7 = 721 and P = 712, respectively. The 
master equation implies an infinite set of equations of motion for expectation values. The 
above-mentioned limitation for the density matrix elements involved, translates into a closed 
system of equations of motion for a reduced number of relevant expectation values. These 
are 



c n = (\l)(l\tf n b n ) , 
v n = {\2)(2\tf n b n ) , 
t n = -ig(\2)(l\tf n b n - 1 ) 



n = 0,1,2... 



n = 0,l,2 



n = 1,2,3... 



and 



(5) 
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Obviously, the average population of the upper (lower) level is given by Co (vo), which 
obey Co + i>o = 1. Of a special interest for the photon statistics are the expectation values 
Vn = (&^ n & n ) = c n +v n , and in particular the average photon number p 1 , which is also denoted 
by N. The imaginary prefactor in the definition of the multi-photon assisted polarization 
t n makes it a real- valued quantity. 

The equation of motion for the expectation value of a given operator A can be obtained 
from the master equation via 



l<^{4" 



llAHjc^ + ^db'.Alb + b^A.b]) 






lii_ 

2 



a\ p A 



°ij + °ij [A CTi 



(6) 



For those in Eq. (jSJ) this leads to the following equations of motion and the corresponding 
steady-state conditions [8|, |l2j, 



121 



{riK + 7)c n + Pv n - 2t n+i = 



d_ 

dt Cn = 

d 

—v n = 7C„ - (nn + P)v n + 2t n+ x + 2nt n = , 
at 

d 2 2 2 (2n-l)K + P + 7 

— t n = g c n + g ncn-i - g v n t r 

at 2 







(7) 
(8) 
(9) 



By adding Eqs. (J7J) and (jHJ) one obtains the steady-state balance relation between the losses 
from the cavity and its feeding through the photon-assisted polarization 



k p n = 2 t n , n > 1 



(10) 



Using this condition and Eq. ((7|) with n = 0, one obtains Pvq = 7C0 + K N, which allows to 
express the steady-state level occupancies in terms of the photon output N 



Co 



P-kN 

P + l 



Vo 



7 + /tiV 
P + l 



(11) 



The unknowns c n and v n can be eliminated from Eqs. (J7J) and © in favor of t n and, 
using again the balance condition Eq. ffTQj) . one is lead[8j to a three-term recursion equation 
for the photonic quantities p n 



A n p n+1 + B n p n -C n p n ^i =0 , n>l 



(12) 



with 



2k 

UK + P + 7 

n k — P + 7 rift (2n — 1) k + P + 7 

"~nK + P + 7 + (n-l)K + P + 7 +K ip ' ^ ' 



(n - 1)k + P + 7 

By using the well-established connection between three-term recursion problems and con- 
tinued fractions [20], Eq. ( Fl2|) allows for obtaining directly steady-state values, without re- 
sorting to the time evolution 17[ . The convergence of the continued- fraction solution is very 
good in all points, except the intermediate pumping region where the transition takes place 
and the population inversion becomes linear. 

In that interval an excellent agreement with the numerical solution can be obtained by 
the following simple ansatz: Assume that (i) Eq. ( TT2|) is valid for n = too, and (ii) the last 
term CoP_i takes in this case the value 0, so that one has 

A oPl + B oPo = 0. (14) 

With p = 1 and the values A and B as in Eq. ffT3|) one obtains for the average photon 

number 

N B P-7 (P + 7 )(P + 7-ftQ n ,, 

N = -A- = —k 8? • (15) 

Using this in Eq. (TTTj) leads for the population inversion w = c — v to 



P + 7- K 



W = K — 13 — . ( 16 ) 



showing that the linear P-dependence of w is a direct consequence of the ansatz. 
In the scaling limit, k — > en, g — > ^Jeg with e — >■ 0, these results become 

P + 7 
W = K ~kf ' (1?) 



for the population inversion and 



N = ^>- ( Z±ll. ( 18) 

2k 8g 2 v ; 



for the rescaled photon output iV = eN. 
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Similar results can be obtained for the three-level model (details of the calculations are 
left for a future publication) with the conclusion that in the scaling limit the population 
inversion behaves in as in Eq. (|T7|) above, while the rescaled photon population obeys 

(P- 7 )y (p + 7 )(p 7 + P y + 7 y) 

k(p + 2V) V(p + 2y) ' l j 

Note that in the limit of large 7 ', that is for very fast |3) — > |1) relaxation, one recovers the 
two-level result, Eq. ( TTBl) . as expected. The agreement of these expressions, Eqs. ( fT7T - [T9l) . 
with the numerical simulations is illustrated in Figs. [2] and |U 

It is clear that the ansatz makes sense only in the interval of P values for which N > 0. 
The first term in Eq. (Tl5|) is positive if P is not too small. On the other hand, for large P 
values, the second, negative term becomes dominant, so that the positivity condition can 
hold only for a finite interval. 

Needless to say, the very good agreement between the ansatz and the numerical results 
(in the interval where the former makes sense) does not constitute a valid proof of the 
former. It is not immediately obvious why Eq. (TT2|) should hold for n = 0. Neither can 



one use C = as an argument to replace the last term C p-i with 0|l9], because it also 
contains the ill-defined p-\. In order to prove the result one has to show first that, indeed, 
the three-term recursion relation can be extended for n = 0, identifying in the process the 
quantity appearing in the role of CqP-i- In a second step, one has to show that, in certain 
conditions, this quantity does vanish, as required by the ansatz. 

The first step of this program is the easier part. It relies on the Glauber-Sudarshan (GS) 



P-representation for the photonic density operator 21] as an integral over the complex plane 
of coherent states 

/j2 
\a)V{a)(a\—. (20) 

Since there is no preferred phase angle in the theory (the reduced photonic density operator, 
obtained by tracing out the emitter indices, is diagonal in the photon number basis) the 
P-function depends only on s = \a\ , V(a) = V(s). Then, the normal-ordered photonic 
expectation values p n turn out to be moments of the quasi-distribution defined by P 

/"OO 

Pn = / s n V(s)ds. (21) 

Jo 

Note, for further reference, that the Poissonian statistics, characterized by p n = X n would 

correspond to a sharp peak in the GS function V(s) = 5(s — A) with A > 0. 

11 



The P-representation contains, in principle, the same information as the density operator, 
but here we take advantage that it allows a natural extension for the definition of expectation 
values. A case in point is p n which, using Eq. ( 121 p . becomes well-defined even for n taking 
continuous, (not just integer) positive values. Using the differential equations which translate 



the master equation Eq. ([2]) into the P- representation formalism J21(, one recovers the three- 
term recursion formula with a continuous index, including the result for index zero 17j|. The 
latter can also be obtained by taking the limit n — > in Eq (fT2|) and using for evaluating 

l^nPn-l 

/»oo /*oo 

np n ^ = / (s n )'V(s)ds = - / s n V'(s)ds ^% P(0). (22) 

Jo Jo 

It is now obvious that the third term in the recursion relation Eq. (|T2|) for n = is propor- 
tional to P(0) and therefore the validity of the ansatz is equivalent to the requirement that 
V(s) vanishes for s = 0. 

The second task is to establish the conditions when the vanishing takes place. This is the 
central point not only in the justification of the ansatz, but also in the proof of the scaling 
limit result. The latter amounts to showing that in this limit the P-function becomes a 
(5-distribution concentrated on a positive value, and this obviously entails that indeed, its 
value at the origin becomes vanishingly small. 

To this end we look at the differential equation obeyed by V(s), and which translates the 
master equation into the langua ge o f the P-representation. We summarize here the main 
steps, the details can be found in 17|]. To start with, the density matrix of our problem has a 
two-by-two block structure, corresponding to the two levels of the emitter. Correspondingly 
one has four P-functions placed in a matrix Vij(s), i,j = 1,2 and the photonic function we 
are interested in is obtained by tracing out the emitter-state indices P = Vx t \ + V%^. Using 



the rules for mapping the master equation for p into a Fokker-Planck equation for P, 21 ] 
one is lead to a system of equations which is the counterpart of Eqs. (jTHH]). After eliminating 
Pi 5 i and p2,2 in favor of P we obtain a second-order differential equation for the latter. 

This equation has then to be analyzed in the scaling limit. To simplify the notation we 
take k itself as the scaling parameter which goes to 0, and impose the condition g — > with 
g 2 /n fixed, by writing g 2 = g 2 n and keeping g 2 constant. The rescaled expectation values 
p n = np n , which are the object of the scaling statement, can be obtained as moments of a 

12 



rescaled P-function 

p n = I K n s n V(s)ds=- I t n V{t)dt, (23) 



OO 



'0 

with t = ks and 

1 _ ft 



nt) = -V[- K ). (24) 

If, indeed, the number of photons increases to infinity in the scaling limit, then the V quasi- 
distribution function moves its weight to larger and larger values and its moments cease to 
exist. In this situation only the rescaled function remains meaningful. Intuitively, according 
to Eq. ff24|) . the graph of the rescaled function V(t) is obtained from that of the original 
V(s) by compressing the latter by a factor of 1/k along the abscissa and expanding it by 
the same factor along the ordinate. This would bring the rescaled function to S(t), in the 
limit k — > 0, were it not for the opposite tendency of V(s) to move away from the origin, 
as discussed above. The net result of these competing trends is what one has to establish. 
It is easy to rewrite the differential equation obeyed by V(s) into the corresponding one for 
V(t), and retain in the coefficients only the dominant terms in the scaling parameter k. The 



result is 



t 2 



1 + p \ 1 
8g 2 J 2 . 



Kp + (t-u)P = 0, (25) 



v= '-l±-iLZ>>-, (26) 



with v an essential parameter in the discussion 

P -l _ (P + 7 ) 2 
2 8~g 2 

whose P dependence is important and therefore sometimes emphasized by the notation 

v(P). Note that v is the same as the rescaled photon population kN, see Eq. (Tl8l) . so that 

they change sign simultaneously. 

The appearance of the small parameter k along with the derivatives suggests a WKB 

approach to the k — V asymptotics of the solution. In other words one searches the solution, 

up to a normalization factor, in the form 

Vit) = exp f~ <p(t)\ , (27) 

in which tp(t) is taken in the leading, zeroth order in k. It is clear that when k gets smaller, 
the value of V(t) around the minimum of <p(t) is greatly enhanced, in comparison with the 
values at other points which, in the view of normalization, become negligible. In the limit 
one obtains a 5-function concentrated at the minimum of (f(t). 
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The equation obeyed by ip(t) in the leading order has the form of a quadratic equation 
for its derivative 



4<^ 2 



7 + P \ 1 
8g 2 J 2 



p' + [t - v{P)\ = . (28) 



Around £ = one of the roots behaves like <p' ~ 2g 2 P/t 2 , i.e. y ~ —2g 2 P/t which in Eq. (j2 



leads to a strongly singular solution. The regular one comes from the other root for which 
yy(0) = -2u(P)/P. 

Two cases arise, depending on the sign of u(P): (i) As long as u(P) is negative, <f'(0) > 0, 
then t = is a minimum for ip(t) and, according to the above discussion, V(t) tends to S(t) 
in the scaling limit, (ii) When v(P) becomes positive, tp' starts at t — with negative 
values and crosses the abscissa at t — v. Then tp(t) has a local maximum at the origin and 
therefore the values of "P(O) become vanishingly small in the limit k — > 0. This is precisely 
the requirement for the ansatz to hold. The function V(t) is now concentrated at the point 
of minimum t = v{P)- 

As a consequence, in the interval in which v(P) is positive one has V(t) — > 8(t — v) and 
all the rescaled expectation values become p n = u n , n > 0. Outside this interval V(t) — > S(t) 
and all p n = 0, except p Q = p Q , which is equal to 1 by definition. The change is abrupt 
and takes place at the interval endpoints defined by the quadratic equation u(P) = 0. The 
condition for this equation to have real roots is g 2 > 27, or g 2 > 2kj, and then the roots 
are both positive 

P ± = 2~g 2 - 1 ± 2gy/g*-2-y. (29) 

The lowest one, P_ = Pthr, corresponds to the onset of lasing and the highest, P + = P sq , to 
self-quenching. The condition g 2 > 27 distinguishes the two behaviors illustrated in FigJT] 
because, when not fulfilled, no transition takes place. With this, the proof of the scaling 
limit is complete. 

It is instructive to see the action of the scaling limit directly on the recursion relation 
Eq. ( 1T2|) . The essential point is the observation that the n-dependence of the coefficients 
A n ,B n , C n gradually disappears in the limit k — >■ 0. More precisely, the recursion for the 
rescaled expectation values 

A 

— pn+l + B n p n - KCnpn^i = , 71 > 1 , (30) 

K 
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FIG. 5. Second, third and fourth correlation function for (a) the two-level emitter and (b) the 
three-level emitter. The parameters are the same as in Figj2]with e = 1 in (a) and e = 0.004 in 
(b). 



reduces, in the k — y limit, to 



Pn+l 



Bo, 

■K-7-Pn= VPn , 

A-0 



n > 1 



(31) 



with the obvious solution p n = v n ~ x f>\. For the values of the pump where v{P) is negative 
only the trivial solution p n = 0, n > 1 is possible, in order to avoid negative results for 
positive expectation values. On the other hand, when v(P) > 0, Eq. (1311) holds for n = too 
and one has p± = vpo = v. Then the solution is exponential p n = z/ n , n > in accordance 
with the Poisson statistics. 

As k approaches 0, the product tlk in the coefficients of the recursion relation vanishes, 
and this is how their n-dependence is lost. In the process it is the low-index coefficients 
that are the first to get close to their limit values, because the limit requires the product 
uk to be small. Therefore the Poissonian condition g^ = 1 is obeyed by g^ first, and by 
g( 3 \g( 4 \ . . . only later. This is numerically confirmed, as seen in Figj5l and shows that using 



9 



(2)> 



1 as a criterion of truly coherent light may be, in this sense, somewhat premature. 



IV. CONCLUSION 

By solving the master equation for a single emitter in JC interaction with a cavity mode 
one observes a sudden change in the behavior of the steady-state solution. This is indicative 
of the onset of lasing and offers the possibility of identifying the conditions for a sharp tran- 



15 



sition to a pure, as opposed to approximate, Poissonian statistics. We have shown that these 
conditions imply an asymptotic regime for the parameters controlling the generation and 
loss of cavity photons. Specifically, the domain of parameters for which a sharp transitions 
occurs is defined by both the cavity loss k and the JC coupling g going to 0, provided that 
g scales like t/k. 

The result is supported by numerical data, as exemplified for a two-level and a three- 
level emitter, and is proven using analytical methods for the two-level model. In a previous 
paper [13], the same scaling limit was shown to give rise to a sharp transition and to reproduce 
the threshold value known in the literature, for the random injection model of Scully and 
Lamb. It should be noted that the Scully-Lamb model does not belong to the class considered 
here: while the latter are "embedded" JC systems the former is rather an "intermittent" JC 
one. This fact, together with the numerical evidence, suggests that our scaling limit result 
has a range of validity that is larger than the set of cases for which a full analytic proof is 
available now. 
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